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Abstract 

We derive and discuss the solution of the Boltzmann equations for leptogenesis in a phe¬ 
nomenologically viable SU(5) x A 5 golden ratio flavour model proposed in 00 The 
model employs, in particular, the seesaw mechanism of neutrino mass generation. We find 
that the results on the baryon asymmetry of the Universe, obtained earlier in using 
approximate analytic expressions for the relevant CP violating asymmetry and efficiency 
factors, are correct, as was expected, up to 20-30%. The phenomenological predictions 
for the low energy neutrino observables, derived using values of the parameters of the 
model for which we reproduce the observed value of the baryon asymmetry, change little 
with respect to those presented in . Among the many predictions of the model we hnd, 
for instance, that the neutrinoless double beta decay effective Majorana mass mee lies 
between 3.3 meV and 14.3 meV. 


^E-mail: julia.gehrlein®student .kit. edu 

^Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, 
Bulgaria. 

^E-mail: martin.spinrath@kit.edu 
■^E-mail: xinyizhangl8@gmail.com 



1 Introduction 

The origin of flavour is one of the most challenging unresolved fundamental problems in 
particle physics. The questions of why there are three generations (not more and not less), 
of the origin of the hierarchies of the fermion masses and of the very different quark and 
neutrino mixing patterns are still far from having received a satisfactory explanation. 

In recent years an approach to the problem of flavour based on discrete flavour symmetries 
became widely used especially in treating the flavour problem in the lepton sector, for a recent 
review see, e.g., [^. A large number of models employing discrete flavour symmetries have 
been proposed. However, many of these models focus only on leptons and only reproduce the 
observed neutrino mixing angles with possibly a few additional predictions for the leptonic 
CP violation phases and/or the absolute neutrino mass scale. 

Here we will focus on a particular model which reproduces all flavour information in the 

quark sector and in addition to reproducing the mixing angles in the neutrino sector, provides 
predictions for the absolute neutrino mass scale and the leptonic CP violation phases. In 
we discussed a slight modification of the original model, which allowed us to accommodate 
successfully the generation of the baryon asymmetry of the Universe within the leptogenesis 
scenario . In that previous publication we used analytic approximations to calculate the 
baryon asymmetry, which can be expected to be correct only up to 20-30%. In the present 
article we go beyond these approximations and calculate the baryon asymmetry by solving 
the relevant system of Boltzmann equations numerically. We show that using this more 
precise method of calculation of the baryon asymmetry one can still generate successfully the 
observed value of it in the model considered. We discuss also the impact of the new results on 
the baryon asymmetry on the predictions of the low energy observables of the model - on the 
correlation between the angles 0 i 3 and 023, on the values of the leptonic CP violation phases, 
on the value of the effective Majorana mass in neutrinoless double beta decay, ruee, etc. 

The article is organised as follows. In Section 2 we review the model constructed in 
and its modihcation proposed in [^ . In Section 3 we discuss the Boltzmann equations and the 
solutions for the baryon asymmetry we obtain. We update the results on the neutrino masses 
and mixing angles previously obtained in in Section 4. Section 5 contains summary and 
conclusions. 

2 The Leptonic Yukawa and Majorana Mass Matrices 

In this section we briefly recapitulate the Yukawa couplings and Majorana mass matrices in 
the lepton sector of the model of interest to fix notations. The structure of these matrices 
is justified by the flavour symmetries of the model and is discussed extensively in [^[^. The 
interested reader is referred to these articles for details. 

The right-handed neutrino Majorana mass matrix reads 




^RR = 2/2 -\/3u2 

y -V^V2 


( 2 . 1 ) 


where V 2 and are complex (vevs) of a flavon breaking the A 5 family symmetry. This matrix 
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is of the golden ratio pattern type A [^, i.e., it is diagonalised by 
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( 2 . 2 ) 


where is the golden ratio and Pq is diagonal matrix containing the two CP violation 

Majorana phases ai and 0 : 2 , Pq = Diag(exp(—^), exp(—^), 1 ). The phases ai/2 and 02/2 
are related to those in the convention used by the Particle Data Group [^, 021/2 and 031 / 2 , 
as follows: 021 = 02 — 01 , 031 = —oi. 

The matrix of charged lepton Yukawa couplings has the form: 


/ 0 -1/2021 0 \ 

Ye = 6 ai 2 6022 6032 , (2.3) 

Vo 0 - 3 / 2033 / 


ij are complex parameters which are fixed by the quark sector (since it is a GUT 


where the a. 

model) and the charged lepton masse s [I] . Note that we did not use here standard GUT 
relations but the relations proposed in [Tpwhich are in good agreement with the current data 
on fermion masses and the Higgs mass results [S-IO . Since the Ojj depend on tan/3 and to 
redo the fit is very time consuming we have fixed this parameter here to 30 which is in good 
agreement with the aforementioned GUT relations. 

The matrix of neutrino Yukawa couplings can be written as 


Y, = yLO + 5 Y, . 

The matrix 

/I 0 0\ 

yLo = 2/H 0 0 1 , 

Vo 1 0/ 

appeared in the original model [^, while 

/O 1 0\ 

6Y, = \y^\cP^ -1 0 0 , 

Vo 00 / 


(2.4) 


(2.5) 


( 2 . 6 ) 


was introduced in the second study |^. The phases can be chosen such that is real and 
then it turned out that the baryon asymmetry is generated only by the small correction 5Yi, 
(note that c <C 1). For more details and analytical estimates the reader is referred to [^. 

It proves convenient for our further discussion to define the parameters which will play an 
important role in the analysis we will perform: 

Mi = ^(Y + Y) = ^|Y||l + pe''^|e''^S 01 = arg(Y + Y) , (2.7) 

M 2 = ^(Y-y) = ^|Y||l-pe'^|e''^^ 02 = arg(Y - Y) , ( 2 . 8 ) 

M 3 = 03 = arg(Y) , (2.9) 
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where 


X = (4^3 + V2)y2 , 

Y = 3'/5v2y2 , 

Y 

X ’ 

(j) = arg(y) — arg(X) . 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


One of the Majorana phases which we choose to be <pi can be set to zero by applying a 
redefinition of the heavy Majorana fields. The remaining two phases (f >2 and (^3 can be 
expressed in terms of p and cj) using the complex mass sum rule Mi + M 2 = M 3 

|M3|2-|Mi| 2-|M2|2 _ l-p 2 

2IM1IIM2I Vl-Vcos2</- + /i 4 ’ 

|Mip - IM 2 P + IM 3 P l + pcoscj) 

cos Cy 3 = -;-r;-;- = - , . 

2 IM 1 IIM 3 I Y'T+^)(7cos^^r)c^ 

We note that only normal ordering is viable in the model considered , and the Yukawa 
couplings are degenerate in LO, so that we have IM 3 I < IM 2 I < |Mi|. The Majorana phases 
«!, 02 are employed in the renormalisation group evolution package REAP |11| we are going 
to use in our numerical analysis, and the phases (/>2 and (/)3 are related, up to corrections of 
order c^, via 


(2.14) 

(2.15) 


oi — —(/>3 and Oi2 — 4’2~4>3- 


(2.16) 


3 Boltzmann Equations 


In this section we discuss the Boltzmann equations for this model. Since we set the lepto- 
genesis scale to be Ms ~ 10^^ GeV, we find with tan/3 = 30, 10®(1 + iav? (3) GeV< Ms < 
10^^ (1 + tan^/3) GeV. Values of Ms in this interval correspond to the two-flavour lepto- 
genesis regime |13[ 14 . We will perform the analysis of the baryon asymmetry generation in 
this regime. 


We use the set of Boltzmann equations in supersymmetric leptogenesis |15 -17 which 


we briefly summarise below (for notational details and further explanations, see the original 
papers). 

The baryon asymmetry generated in the two-flavour regime in leptogenesis is determined, 
in particular, by the evolution of the heavy Majorana neutrino and sneutrino number densities 
(abundances), and , and of the lepton charge and CP violating asymmetries in the 

charges Lg + and Lr, Yas = + Ya^ and Y/s,^ (where Ya^ = Ya^ -|- Ya^-, with = 

B/3 — during the epoch of the evolution of the Universe when the abundances Y^r. and 

IV start to deviate from their equilibrium values (the out-of-equilibrium Sakharov condition 
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The evolution of the quantities of interest in the epoch of interest can be described 
by a system of coupled Boltzmann equations, which in the case of the MSSM and of the 
two-flavour regime we consider read 


15-17 
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(3.4) 


In these expressions, Y^^'^ is the Ni (equilibrium) abundance, 2 : = Myi,g/T with Mflg 
being the mass scale of flavoured leptogenesis and T being the temperature of the thermal 
bath, s = (7*27r^T^/45 is the entropy density, H(T) ~ 1.66y^r^/mpi is the expansion rate 
of the Universe, g* = 228.75 and the Planck mass mpi ~ 1.22 x 10^® GeV. 7 }^ is the thermally 
averaged Ni decay rate, 7 ^^=! is the AL = 1 scattering rate of Ni with leptons, quarks 

and gauge bosons, flavour dependent N inverse decay rate, and 7 ^al=i 

washout rate of the AL = 1 scatterings. The A matrix that relates the B/3 — Li asymmetry 
1 a; and the lepton charge asymmetry in doublets, Y/, defined in Y; = A/pYAj,, in MSSM 
and the two-flavour regime reads 


1 /-541 152 \ 

761 V 46 -494y 


(3.5) 


The expressions for the CP violating asymmetries generated in the heavy Majorana (s) neutrino 
decays are the same as those derived in [^: 
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Figure 1: Solutions for the Boltzmann equations using the example input set. The horizontal 
grey band represents the 3cr region for the observed baryon asymmetry Yb = (8.58 ± 0.22) x 
10“^^, where, for simplicity, we obtained the Scr region by multiplying the Icr error by a factor 
of three. The dot-dashed line is See text for further details. 
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Due to the fact that leading order neutrino Yukawa coupling is unitary (except for an overall 
factor (y”)^), we have e? = ef + = —ej to leading order. 

We do not include here the thermal corrections to the CP asymmetries since they give 
only negligible contributions within the low temperature regime 19,20 . We neglect also the 
AL = 2 process for the same reason as stated in j^, i.e., given the value of y”, the AL = 2 
processes do not have a significant impact at the leptogenesis scale of interest Ms — 10 ^^ 
GeV. 

The final baryon asymmetry is 


Yb=^^{Ya,+Ya^). (3.9) 

As an example we show the result of the Boltzmann equations in Fig. for a single set of 
parameters chosen from the numerical scan in the next section. The parameters for this plot 
are 


tan/3 = 30, Mg = 8.51 x 10^^ GeV, c = 0.053, p = 6.29, 7 = 3.33 . (3.10) 
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Table 1: The best-fit values and the Scr ranges for the parameters taken from 
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The two 


minima for both 613 and 023 correspond to normal and inverted mass ordering, respectively. 


And we have 


Yb = 8.86 X 10"^^ , (3.11) 

= -4.72 X 10“^° , (3.12) 

y^^ = 7.47 X 10“^° . (3.13) 

Note that at this point we had given two different example points in our previous publi¬ 
cation [^. But with the improved calculations none of them is in good agreement with the 
experimental result on Yb anymore and hence we have chosen here a different example point. 
This shows furthermore the importance of the current study. 

4 Results 

In this section we show the results for the masses and the mixing angles as well as the results 
for Yb in dependence on the parameters in our model. To obtain realistic values for Yb 
coming from the numerical solution of the Boltzmann equations given in Sec. we had to 
slightly increase the value of the neutrino Yukawa coupling from yi = 0.1 to yi = 0.12 which 
is a 20% change as anticipated. 

4.1 Masses and mixing angles 

For our numerical scan we use the method as described in [^. For the parameters which 
characterise the charged lepton and quark sector we use the fit results given there with 
tan (3 = 30 and Msusy = 1 TeV. The renormalisation group evolution is done using the 
REAP package [IT] . 

There are six free parameters left in our model, the moduli \X\ and |Y|, the phases (p 
and b \2 in the leading order matrices and the modulus |c| and the phase 7 , which come 
from the correction to the neutrino Yukawa matrix. We performed a random scan over 
these parameters and impose the experimental ranges for the mixing angles, mass squared 
differences (cf. Tab. and Yb as constraints. For Yb we used 

Yb = ( 8.58 ± 0.22) x 10“^^ , (4.1) 


6 







where the 3a uncertainty is obtained, for simplicity, by multiplying the la error by a factor of 
three. In order to calculate Yb we solved the Boltzmann equations given in Sec. [^numerically. 

Before we present the results for the normal ordering of the neutrino masses we comment 
briefly on the case of inverted ordering. In the original model the inverted ordering was 
not viable due to incompatible constraints for 9 i2 coming from the mass sum rule on one 
hand and from the angle sum rule in our model on the other hand. Due to the correction for 
the neutrino Yukawa matrix the mixing angles are modified with corrections of order c which 
nevertheless have to be of the order of c ~ 0.4 to save the inverted ordering as we have shown 
in an estimate given in [^ . Since c is associated with a higher dimensional operator such high 
values are not plausible. 

Turning now to the normal ordering, we show the results for the masses and mixing angles 
in Fig. [^ We find all parameters in agreement with the experimental la (3cj) ranges for the 
masses, mixing angles and Yb- 

The correlations between @13 and the Majorana phases are weaker than in [^ which is due 
to the constraints on the phase (p and the value of \Y\ coming from the numerical solution 
of Yb- Furthermore we need here a larger value of c which washes the correlations out and 
enlarges the ranges for the phases. Nevertheless we find the phases to be in similar ranges as 
in [^, namely we obtain 


<5 G [9°, 119°] or [239°, 344°] 
ai G [0°,134°] or [220°, 360°] 
02 G [66°, 134°] or [222°, 282° 


(4.2) 

(4.3) 

(4.4) 


For the rephasing invariant Jcp which determines the magnitude of CP violation effects in 
neutrino oscillations 26 we find values in the ranges Jcp = ±(0.006,0.036). Our predictions 
for nentrinoless double beta decay are shown in Fig. [^ For the lightest neutrino mass mi 
which is mostly determined by the mass sum rule we obtain values between 12 meV and 
22 meV. For the observable in neutrinoless double beta decay niee we obtain values between 
3.3 meV and 14.3 meV. For the sum of the neutrino masses we predict 

^ m^ G 0.077 - 0.099 eV , (4.5) 

which might be determined, e.g., from cosmology. So far there is only an upper bound [ 22 ] 

^m^<0.23eV, (4.6) 

which is well in agreement with our prediction. The second observable is the kinematic mass 
mjs as measured in the KATRIN experiment 27 which is given as 


2 222 , 222 , 22 

m^ — ™'lCi2Ci3 ± 1 ^ 12 ^ 12^13 + "^3'®13 


(4.7) 


Here we predict m^ ~ 0.015 — 0.023 eV which is below the projected reach of m^ > 0.2 eV 
of KATRIN. 


4.2 Leptogenesis 

Moving on to our predictions for Yb, we show the results of our parameter scan in Figs. 
and [H 
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Figure 2: Results of our numerical parameter scan. Blue (red) points are in agreement within 
3cr (Icr) of the low energy neutrino masses and mixings and Yb in our model. The allowed 
experimental 3 ct (Ict) regions are limited by blue (red) dashed lines. The black dashed lines 
represent the Icr range for the not directly measured CP phase <5 from the global fit 
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Figure 3: Prediction for the effective neutrino mass niee accessible in neutrinoless double beta 
decay experiments as a function of the lightest neutrino mass rrii. The allowed experimental 
3 ct (Itr) regions for the masses and mixing angles in the case of normal ordering are limited by 
blue (red) dashed lines. Blue (red) points are in agreement within Scr (Ict) of the low energy 
neutrino masses and mixings and Yb in our model. The grey region on the right side shows the 
bounds on the lightest mass from cosmology 22 and the grey region in the upper part displays 
the upper bound on the effective mass from the EXO experiment 24 . The red, straight lines 


represent the sensitivity of GERDA phase I respectively GERDA phase II 25 
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Figure 4: Results of our numerical scan for the total baryon asymmetry Yg in dependence 
of the six most relevant parameters. Blue (red) points are in agreement within 3a (Icr) of the 
low energy neutrino masses and mixings and Yb in our model. 
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Figure 5: Results of our numerical scan showing the correlations between the lightest right- 
handed neutrino mass and the six most relevant parameters for leptogenesis. Blue (red) points 
are in agreement within 3(t (Ict) of the low energy neutrino masses and mixings and Yg in our 
model. 
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We find parameter points which are in agreement with the la (3a) range for Yb- The 
ranges for the (p and p are as in our previous study between 4.7 and 7.6 for p and between 
0.66 and 1.38 or -1.38 and -0.66 for p since they are mostly determined by the masses and 
mixing angles. However the constraints coming from the numerical solution of Yb now forbid 
values of p between -0.66 and 0.66. Interesting is also the range for the mass of the lightest 
right-handed neutrino which is in the range from 0.778 • 10^^ GeV to 0.862 • 10^^ GeV which 
shows a clear correlation to c but no correlation to 7 .The other correlations can be seen in 
Fig. [51 We furthermore want to note that the rather narrow range of M 3 is due to the fact 
that we fixed y” to a certain value. As a rule of thumb the scale of light neutrino masses has 
to remain the same and hence a variation of 10 % in y” leads to a variation of 20 % in M 3 . 

The size of the correction c has to be in the range between 0.035 and 0.085 and 7 is in 
narrow ranges around 0 or tt. This also explains the disconnected regions for the mixing 
angles and phases which did not appear so in our previous studies. As in our previous study 
the sign of Yb is uniquely determined by 7 . For the Icr ranges of the parameters the values of 
7 have to be in the region around 0. This implies that we only find parameter points for the 
Icr ranges compatible with a negative value of Yg. Note that this is mainly driven by the Icr 
range of 823 which is not very well determined. Indeed, with the results from the Valencian 
fitting collaboration 28 where 623 is allowed to be in the second quadrant, even at the la 


level positive values for Yb are possible. 

As 7 is in narrow ranges around the CP conserving values, the Majorana phases are the 
major sources for CP violation here, as it can be seen from Fig. [^ 


5 Summary and Conclusions 

In this work we have revised predictions for leptogenesis in an SU(5) x A 5 golden ratio GUT 
flavour model from [^. In that publication we had used approximations to calculate the 
baryon asymmetry which are known to be precise only up to 20-30%. Instead we have 
solved here the full set of Boltzmann equations numerically and could show that we can still 
successfully accommodate the experimental values of all mixing angles, fermion masses and 
the baryon asymmetry even at the la level. Nevertheless, to do so we had to adjust some 
parameters. For instance, we have increased the neutrino Yukawa coupling by 20% from 
yi = 0.1 to yi = 0 . 12 . 

All of the main features of the original model are still valid. For instance, the neutrino 
mass sum rule is correct up to small corrections and the inverted ordering is still ruled out. 
Furthermore, the Yukawa coupling ratios Hr/yb ~ —3/2 are unaffected by the modification 
of our model and hence imply non-trivial constraints on the spectrum of the supersymmetric 
partners of the Standard Model particles. 

The phenomenological predictions are not changed significantly as well. For the lightest 
neutrino mass we And the range mi G [0.012,0.022] eV, but the strong correlation between O 13 
and the phases is somewhat weakened due to an increased value of c. For the phases we find 
5 G [9°, 119°] or [239°, 344°], ai G [0°, 134°] or [220°, 360°] and 02 G [66°, 134°] or [222°, 282°]. 
Related to the Majorana phases is the observable in neutrinoless double beta decay where we 
predict m^e between 3.3 meV and 14.3 meV. 

We have proposed the first SU(5) x A 5 model to our knowledge which can successfully 
accommodate for all parameters in the neutrino sector and for the baryon asymmetry of the 
Universe as calculated from the Boltzmann equations and for the not yet measured quantities 
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in the neutrino sector we make testable predictions. For this our model needs comparatively 
few parameters which makes it very appealing and testable in future experiments. 
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